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ABSTRACT 

We present an algorithm (MEAD, for 'Mapping Extinction Against Distance') which will de- 
termine intrinsic (/ — i') colour, extinction, and distance for early-A to K4 stars extracted from 
the IPHAS r' /i' /Yia photometric database. These data can be binned up to map extinction in 
three dimensions across the northern Galactic Plane. The large size of the IPHAS database 
(~ 200 million unique objects), the accuracy of the digital photometry it contains and its faint 
limiting magnitude (/ ~ 20) allow extinction to be mapped with fine angular (~ 10 arcmin) 
and distance (~ 0.1 kpc) resolution to distances of up to 10 kpc, outside the Solar Circle. 
High reddening within the Solar Circle on occasion brings this range down to ~ 2 kpc. The 
resolution achieved, both in angle and depth, greatly exceeds that of previous empirical 3D 
extinction maps, enabling the structure of the Galactic Plane to be studied in increased detail. 
MEAD accounts for the effect of the survey magnitude limits, photometric errors, unresolved 
ISM substructure, and binarity. The impact of metallicity variations, within the range typical 
of the Galactic disc is small. The accuracy and reliability of MEAD are tested through the 
use of simulated photometry created with Monte-Carlo sampling techniques. The success of 
this algorithm is demonstrated on a selection of fields and the results are compared to the 
literature. 

Key words: surveys - methods: miscellaneous - ISM: dust, extinction - ISM: structure - 
Galaxy: disc - stars: general 



1 INTRODUCTION 

The INT/WFC Photometric Ha Survey of the northern Galactic 
Plane (IPHAS; iDrew et aLll2005h is the first comprehensive digital 
survey of the disc of the Galaxy {\b\ ^ 5°), north of the celestial 
equator. Imaging is performed in the i' and Ha bands down to 
/ ~ 20 (10a). The high quality photometry from this survey has 
been shown to enable early-A type stars to be select ed very easily, 
due to their strong Ha absorption (iDrew et al.ll2008h . In this same 
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study, the (/ — i') colours provided by IPHAS were used to deter- 
mine the extinctions of the A stars selected, and then distances were 
deduced using photometric parallax. This capability is a specific 
example of a more general property of the IPHAS: since {r' — Ha) 
correlates strongly with intrinsic (/ — i') colour, the combination of 
i' and Ha can be used to disentangle intrinsic colour and extinc- 
tion unambiguously. The aim of this paper is to present a method 
which extends the approach of Drew et al.j (l2008l) to incorporate a 
broader range of spectral types and, therefore, many more objects. 
Then extinction can be successfully mapped at fine angular resolu- 
tions out to large distances. 
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The IPHAS catalogue contains high quality photometry for 
~ 200 million objects, equating to over 10^ objects per square de- 
gree. Our analysis of a wide variety of IPHAS fields has shown that 
the data quality allows for accurate determinations of distances and 
extinctions to the studied objects, with median errors of typically 
~ lOOpc on distance and approximately 0.1 magnitudes of error on 
the determined extinction. The limiting magnitude of the IPHAS 
catalogue (/ ~ 20, 10a) is sufficient for extinctions to be measured 
to di stances up t o ^10 kpc from the Sun. 

Juric et al ] (l2008h discuss and demonstrate the potential for 
accurate wide area surveys in analysing the structure of the Milky 
Way. They utilise the Sloan Digital Sky Survey (SDSS) in order to 
determine the stellar number density distribution mostly at Galac- 
tic latitudes |^| > 25°. As a result of this they are able to eas- 
ily study the Galactic halo and local scale heights of the Milky 
Way's thick and thin discs, and find localised over-densities. In 
contrast the IPHAS survey area covers low Galactic latitudes, en- 
abling the Galactic thin disc to be studied in detail, well beyond 
the solar neighbourhood. The cost associated with this is that ex- 
tinction, which is largely negligible in the SDSS survey area, be- 
comes a significant concern to the extent that understanding its 
three-dimensional (3D) distribution is a necessary hurdle that must 
be overcome on the path to studying the structure of the Galactic 
disc. 

Applications of extinction mapping go beyond studies of 
Galactic structure: a distance-extinction relationship, combined 
with an estimate of the extinction to an object has often been ex- 
ploited to give an estimate of the distance to individual objects. 
Such estimates are of particular use for objects such as planetary 
nebulae, for which reliable estimates of distance are difficult to ob- 
tain (Viironen et al. in preparation). 

There is a long history associated with trying to map the dis- 
tribution of extinction, soins back as far as iTrumplej (Il930h and 
Ivan de Kamd hQSO") who both appreciated that extinction is con- 
centrated near the Galactic plane. Subsequent efforts to map ex- 
tinction can be broken down into two broad categories: 2D maps 
of the asymptotic extinction and 3D maps with varying degrees of 
model input. 

Several studies have produc ed 2D maps of the ma ximum 
asymptotic extinction over the sky. iBur ste in & HeilesI il982) used 
HI mapping and galaxy counts, es timating t he errors in derived ex- 
tinctions to be ^ 10%. ISchlegeL Finkbeiner & DavisI (tl998) pro- 
vided updated maps based o n COBE and FIR AS measurements of 
dust emission at lOO^m. The' Schlegel etal.U l998') maps have since 
been widely used, but they warn that their values for \b\ < 5° are 
su bject to significant e rrors, an as sertion that has been supported 
by lArce & GoodmanI (Il999h and ICam bresv, J arrett & Beichman 
(l2005h . More recentlv lproebrich et al.l 12005) have used cumulative 
star counts from the 2MASS survey to construct a 2D extinction 
map of the Galactic plane (\b\ < 20°). 

3D extinction models have grown progressiv ely more sophis- 
ticated, moving from the simple cosecant model ( Parenago "1945*) 
to the modern, high resolution models of Mendez & van Altena 
([T998h: IChenet al l (Il999l): iDrimmel & Spergell (1200 ll) and 
lAmores & Lepin J "(l2005h . [Marshall et al.l tom compared the 
Besangon Galactic model (iRobin et al.ll2003h . without extinction, 
to 2MASS (Two Micron All-Sky Survey) observations. The dif- 
ferences between the observed and theoretical colours of predomi- 
nantly K & M giants were then attributed to extinction. This semi- 
empirical method is able to map extinction at fine angular resolu- 
tions (0.25° X 0.25°). 

Historically, attempts to empirically map extinction in three 



dimensions have been constrained by the availa bility and qua! 
ity of photo n ietric data. Early studies include 
'Eggen (1963); "isserstedt & Schmidt-Kaler| (Il964h : 
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Neckell 
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Scheffler (1966, 1967), all of which reached no fainter than V ~ 10 
and therefore could only cover relatively limited ranges of up to 
a 1 — 2 kpc. Also, they all contain fewer than 5000 objects, limit- 
ing th eir mappings to very coarse spatial resolutions. Furthermore, 
Fitzgerald (1968) found that these works produced 'varying, and 
somet i nies apparently con flicting results'. Subsequently, Fitzgerald 
(ll968h : lL^ck3 (Il978n and lNeckel. Klare & Sarcandej ( ll98Qh con- 
tinued, with similar approaches, but they were still limited by the 
available pho tometric catalogues, with the most recent of the three 
(iNeckel etal.l l 1980) using a catalogue of ~ 11,000 objects. As a 
result, they achieved angular resolutions on the order of a few de- 
grees. 

Here, we present a new approach to mapping the 3D extinction 
within the plane of our Galaxy, based on IPHAS photometry. The 
depth, quantity and quality of this data set allows us to construct a 
new extinction map with high angular (10^ and distance (100 pc) 
resolution out to very large distances (up to 10 kpc). In the present 
paper, we focus on the algorithm used to construct this map, leaving 
the map itself to be presented in a separate publication. The algo- 
rithm which is used to calculate extinction distance-relationships 
from IPHAS data (hereafter referred to as MEAD, Mapper of Ex- 
tinction Against Distance) is discussed in detail in Section [3] In 
order for it to be used with IPHAS data it is necessary to calibrate 
this translation of IPHAS photometry into derived quantities (Sec- 
tion |2]). Section]?] shows the use of simulated photometry to deter- 
mine the values of some key parameters used in MEAD and demon- 
strate its accuracy in the face of a range of possible problems. Fi- 
nally, Section [5l investigates the use of MEAD on real IPHAS data, 
showing the power of its approach. 



2 DECODING INTRINSIC COLOUR, LUMINOSITY 
CLASS AND EXTINCTION FROM IPHAS 
PHOTOMETRY 

In order to exploit IPHAS, it is necessary to understand the be- 
haviour of its colour-colour and colour-magnitude planes - specif- 
ically the intrinsic col ours of normal st ars and their response to 
extinction. To this end lDrew et al.l (l2005h produced simulations of 
the luminosity class sequenc es on th e IPHAS colour-colour plane 
using the spectral library of IPicklesI (T^98). In Fig. [T] we provide 
an example of (/ — Ha, / — i') data for a survey field with differ- 
ently reddened simulated main sequences spanning the 09 to M3 
spectral type range superposed: this shows how reddening shifts 
the main sequence across to the right, creating a lower edge to the 
main stellar locus that is well-described as the early-A reddening 
line. Reddening also moves O and B stars to the same colours as A 
and F stars, though they are still separated in apparent magnitude. 
For the purpose of emulating the colour behaviour of normal stars, 
the Pickles (1998) library contains limited numbers of objects that 
cannot support a full exploration of the available parameter space. 
In particular, in order to examine the impact of changing metal- 
licity and fine variations in luminosity class, we have recomputed 
expected IPHAS colours from grids of model atmospheres. 

New luminosity class sequences have been obtained for 
the IPHAS colour-c olour and colour-magnitude planes usin g the 
iMunari et alj (l2005l) library of ATLAS9 synthetic spectra. iKuruczl 
( 2005) discusses some particular limitations of ATLAS9 models, 
some of which are relevant to this study: ATLAS9 spectra include 
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Figure 1. Black points show data from IPHAS field 4199, where 13 ^ / < 
20. Main sequences where extinctions equivalent to Ay = (red), 4 (green), 
8 (blue) for an AOV star have been applied are shown with solid lines and 
the dashed cyan line shows an A3V reddening line. There are two main 
populations visible in Fig.[T] the first lies at bluer (r' — i') aligned with and 
extending beyond the simulated unreddened main sequence: this population 
is mainly K-M dwarfs. The main body of stars falls between the Ay = 4 and 
Ay = 8 main sequences. The long trail of objects stretching from {r' — i')^ 2 
are M giants. Finally, there are five objects that appear to lie above the 
unreddened main sequence, these objects are candidate Ha emitters. 
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Figure 2. A Hertzsprung-Russell diagram of the simulated luminosity 
classes used in this study. The unshaded region of the diagram indicates 
the intrinsic colour range of interest, which is the range of intrinsic colour 
that MEAD works on. 



only an incomplete sample of diatomic molecules and, with the ex- 
ception of H2O, no triatomic or larger molecules. This inherently 
leads to deviations fro m reality, which beco me significant for late 
K and all M type stars (iBertone et al.l l2004Y In particular the ab- 
sence of VO lines means that Kurucz models are not available for 
Tgff < 3500^ (approximately M4 or later) as they would be particu- 
larly inaccurate ( K urucz^.200 5). Similar problems affect other grids 
of atmospheres, with the consequence that, for our purposes, there 
are no reliable alternative model grids at these low temperatures. 



Table 1. The intrinsic {r' — i') and (r' — Ha) colours, and r' magnitude com- 
pared to approximate spectral type an d Teff for s olar m etallicity dwarfs. The 
data are derived from the Straizvs & Kuriliene ("1981') reff, l ogj^ calibration 
and the absolute magnitude calibration of Houk et al^ ( Il997h . 



Intrinsic Intrinsic r' Absolute Spectral Teff 



(/ - i') 


(r — Ha) 


Magnitude 


Type 




-0.047 


0.029 


0.12 


B8 


11510 


-0.028 


0.013 


0.51 


B9 


10400 


-0.006 


-0.002 


0.80 


AO 


9600 


0.001 


-0.005 


1.09 


Al 


9400 


0.013 


-0.008 


1.28 


A2 


9150 


0.025 


-0.008 


1.47 


A3 


8900 


0.060 


0.006 


1.84 


A5 


8400 


0.097 


0.032 


2.21 


A7 


8000 


0.171 


0.098 


2.64 


FO 


7300 


0.203 


0.128 


2.92 


F2 


7000 


0.258 


0.174 


3.37 


F5 


6500 


0.296 


0.202 


3.83 


F8 


6150 


0.318 


0.218 


4.11 


GO 


5950 


0.336 


0.229 


4.39 


G2 


5800 


0.372 


0.252 


4.75 


G5 


5500 


0.406 


0.269 


5.21 


G8 


5250 


0.436 


0.283 


5.58 


KO 


5050 


0.454 


0.291 


5.76 


Kl 


4950 


0.472 


0.298 


6.04 


K2 


4850 


0.504 


0.312 


6.20 


K3 


4700 


0.527 


0.322 


6.48 


K4 


4600 



2.1 Method and underpinning stellar data 

So as to use the synthetic spectral libraries to create lu- 
minosity class sequences, the classes of the MK system 
(iMorgan. Keenan & Kellmaniri943h must be linked to Tgff, logg 
and absolute magnitudes {My). This study adopts the Tgff and log^ 
calibrations of IStraizvs & Kuriliend (Il98lh . as well as their abso- 
lute magnitude ca librations fo r all cl asses except dwarfs, for which 
we use the data of iHouk et all (Il997h based on HIPPARCOS paral- 
laxes. In Table [1] we give the solar-metallicity mapping of intrinsic 
{r' — i') colour onto main sequence absolute magnitude and spec- 
tral type that we have adopted, while Fig. [21 shows the relevant 
Hertzsprung-Russell diagram of the five luminosity classes we use. 

All synthet ic colours were determined largely as by 
IPrew et aTl (l2005h . with the modest difference that the atmospheric 
transmission curve at the INT is included in the calculation, so 
that the total transmission profile is the product of the atmospheric 
transmission, CCD response and filter transmission curves. 

As IPHAS magnitudes are based on the Vega system, an 
appropriate spectrum of Vega was required to provide the zero- 
point flux. So as to match spectral resolution (which is neces- 
sary to produce accurate Ha magnitudes) the ATLAS9 Vega model 
spectrum available via kurucz.harvard.edu was used. This assumes 
Teff = 95507^ and logg = 3.80. It should be noted that there has 
been much debate over wh ich ATLAS 9 spectra best fits observa- 
tions of Vega. For example, ICohen et all Il992) and Bohlin (2003) 
argue that the 9400K model (as opposed to the 95 5 OK model used 
in this work) is better suited. However, the differences in derived 
colours are small, 5(/ — Ha)?^ 0.003, so the choice from these two 
models is of little consequence. 

For comparison to the synthetic sequences derived in this 
study, sequences based on libraries of intrinsic (i.e unreddened 
or dereddened) empirical SEDs are employed. At the time of 
writing there are over 100 electronic libraries available, sum- 
marised in the Asiago Database of Spectrophotometric Databases 
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Figure 3. The simulated main (solid black line) and giant (dashed green 
hne ) sequen c es com pared to solar r netallicity main sequenc e objects from 
the IPicklesI (Il998h and STELIB jLe Borgne et"ZI l20Q3h libraries (red 
crosses). As with figure El the unshaded region of the diagram indicates 
the intrinsic colour range of interest, used by MEAD. 
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Figure 4. ZAMS from the IPietrinferni et al.l (l2004h library on the colour- 
colour plane, where the colours have been determined from the appropriate 
iMunari et al] ( l2005i) models. In black is a ZAMS for [Fe/H]= 0, in red a 
ZAMS for [Fe/H]= -1.00 and in green a ZAMS for [Fe/H]=-1.82. 



(ISordo&Munarill200d http://web.pd.astro.it/adsd). However, of 
these only a small proportion cover the wavelength range required 
and many of those are unsuitable. Once these have been filtered 
out, the STEL IB library (iLe Borgne et alJl2003h and the library of 
|Picklesl(ll998h remain. 



2.2 Comparisons between colour sequences from synthetic 
and empirical libraries 

Fig. [3] co mpares the M a in Seq uence derived from the synthetic 
spectra of iMunari et iD (l2005h . to colours derived from empiri- 
cal spectra. For most spectral types the results from empirical and 
synthetic spectra agree well. However, for M-type stars there is 
some disagreement. This could be a result of difficulties associ- 
ated with modelling late type stars, as noted earlier. A compari- 
son was also performed to sequences derived from t he MARCS 
(iGustafsson et al. 2003) and PHOENIX (iBrott & Hausch ildt 2005) 
codes. The MARCS library only covers spectral types from early- 
A to early- K, but it is in excell ent agreement with the results ob- 
ta ined from Munari et al.l (2005b . On the other hand, the sequences 
of 'Brott & Hauschildt ( 2005) are significantly different for spectral 
types later than roughly K4. On further examination the differences 
are almost entirely limited to behav iour in the Ha band. It is worth 
noting that lGustafsson et aTl (l2008h demonstrate that a new version 
of the PHOENIX code, employing the Kurucz line data, performs 
considerably better than the version employed here. Therefore it is 
possible that the discrepancies we encounter will not be present in 

this newer version 

iBertone et al.] (l2004l) finds the ATLAS9 and PHOENIX mod- 
els which best fit intrinsic spectra of known type and list t he T^ff of 
the fit. By assuming that the Tgff calibration of lStraizvs & Kuriliend 
(Il98lh represents the mean temperature for any given class, it is 
possible to calculate the standard deviation of the residuals of the 
Bertone et al. I (l2004h estimates. This is then an estimate of the com- 
bination of the error on the Straizys & Kuriliene (1981) Tgff calibra- 
tion and the natural width of each class. However, allowing for this 
temperature variation makes no measurable difference to our model 
sequences. 

Unfortunately iBertone et al.l (|2004|) do not give the values of 



log g for each fit, so it is not p ossible t o repeat this procedure for 
surface gravity. However. iMunari et al ] (12005) suggest that the nat- 
ural width of each class in surface gravity is ~ 0.25 dex. Applying 
this spread results in the luminosity class sequences changing from 
lines to bands in colour-colour and c olour -magnitude space. The 
majority of objects from the Pi cklesI d 19981) and STELIB libraries 
now fall within these bands. . Bertone et al.l 020041) also demonstrate 
that the quality of the best fit to the stars in the empirical libraries it 
uses, decreases for later types. This behaviour is again attributable 
to the poorly understood molecular lines, which dominate the spec- 
tra of late- type stars. 



2.3 The impact of rotational broadening 

As might be expected, altering the rotational velocity assumed in 
the models has negligible effect on the colours. This is as varying 
rotational velocity will only alter the shape of a non- saturated line 
and will not affect the strength of absorption or emission. Because 
all the filters used are broad compared to emission or absorption 
lines, simply altering the line shapes does not change the observed 
colours. In saturated lines the equivalent width is affected by rota- 
tional broadening, but Ha is the only atomic line which produces a 
significant impact on the photometric colours and it does not satu- 
rate, even in A3V stars where absorption is strongest. 



2.4 The impact of metallicity variation 

Given that IStraizvs & Kuriliend (Il98lh do not consider the effect 
of metallicity it is clearly not possible to use the sequences based 
on it to study the effect of metallicity. However, the library of 
IPietrinferni et al.l (l2004h does consider it and so is used in this sec- 
tion. For ^Munari et all (l2005h models altering metallicity has little 
effect on IPHAS colour-colour diagrams, a shown by Fig.|4] 

Metallicity can, however, have a significant effect on the ab- 
solute magnitude of all objects. Fig. [5] demonstrates this for a zero- 
age main sequence (ZAMS). As it is not possible to determine the 
metallicity of the object from IPHAS observations, it is therefore 
also not possible to uniquely determine the absolute magnitude of 
an object. 
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Figure 5. ZAMS from the IPietrinferni et alJ l2004h library on the colour- 
magnitud e plane, wher e the c olours have been determined from the ap- 
propriate iMunari et S] ( l2005h models. The top, in black, is a ZAMS 
for [Fe/H]= 0. The curve is lowered with decreasing metallicity, in blue 
is [Fe/H]= —0.40, red [Fe/H]= —1.00 and at the bottom, in green, is 
[Fe/H]=-1.82. 
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Figure 6. Main sequences where extinctions equivalent to Ay = (black, 
left), 4 (red, middle), 8 (green, right) for an AOV star have been apphed. 
The dashed black lines show the loci of A3V (bottom), G5V (middle) and 
M4V (top) stars under increasing extinction. 



iFriel et aP ( l2002h measured metallicities in low Galactic lati- 
tude open clusters, fi nding no objects with [Fe/H]< —0.94, whilst 
iCarraro et al determine the mean metallicity in the outer 

disc to be [Fe/H]:^ —0.35. So, observing an object with a metallic- 
ity below [Fe/H]= —0.5 in IPHAS would be unusual and metallic- 
ities as low as [Fe/H]= —1.82 would be extremely rare. Therefore, 
the effect of metallicity variations on absolute magnitude will be 
relatively small in IPHAS observations. Section l48l further exam- 
ines the impact of undiagnosed metallicity variations. 

2.5 Extinction & reddening 

It is very important to note that, in the IPHAS colour-colour plane, 
the reddening vector is at a significant angle with respect to the 
majority of the main sequence (see Fig. [6]). This makes it possible 
to derive accurate estimates of intrinsic (/ — i') colour and extinc- 
tion simultaneously. For colours based on other commonly used 
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Figure 7. The c o efficie nt bj.' against initial (/ — i') for the 
IStraizvs & Kurilien j (Il98ll) main sequence. Note that the coefficients 
a/ and c/, multiplying the squared and constant terms, are roughly 10^ 
times smaller than b/ . 



filter sets such as u' g'r'i'z! or JHK, the reddening vector is almost 
parallel to the main sequence, creating degeneracy in the available 
(intrinsic colour, extinction) solutions. This property of the IPHAS 
colour-colour plane is a result of the strong intrinsic (/ — i') colour 
sensitivity of the {r' — Ha) colour, whilst its response to reddening 
is weak compared to that of broad-band colours. The former prop- 
erty is due to (r' — Ha) being a proxy for the Ha equivalent width, 
which itself is strongly dependent on intrinsic (/ — i') colour. 

In many applications in astrophysics it is acceptable to treat 
the ratio Ry = Ay /E{B — V) as constant across a wide range of 
objects, so that extinction can easily be calc ulated fr om reddening 
and vice-versa. However, as pointed out bv lMcCalll /2004) this is 
not in general true, as the extinction integrated across a given broad- 
band filter is a function of the SED of the observed object and the 
amount and wavelength dependence of the extinction it has expe- 
rienced. iMcCalll (I2OO4I) suggests that monochromatic measures of 
extinction should be used, in preference to broadband measures, as 
they are independent of the observed objects' SEDs. Following this 
concept, MEAD measures a monochromatic extinction at 6250A 
(^625o)- This is near the effective wavelength of the / filter in the 
INT system (i.e. taking into account the WFC and atmospheric re- 
sponse). If the intrinsic colour of the observed object is known it 
will then be possible to convert to and from more commonly-used 
broadband measures. In this study, we will refer to monochromatic 
extinctions in discussions of the simulated photometry, while ex- 
tinctions derived from real data that are compared with the litera- 
ture will be presented as broadband values, A^/. The extinction at r' 
is preferred to the extinction in the V band for the reasons discussed 
bv lCardelli. Clayton & Mathisl(ll989h . Section l4l0l investigates the 
ability of MEAD to change between monochromatic and broadband 
measures of extinction as R is varied. By choosing to normalise the 
extinction law to a wavelength in the middle of the range studied 
the systematic error produced by uncertainty in the value of R is 
minimised. 

This work adopts the extinction laws of 'Fitzpatrickl (l2004l) for 
the calculation of extinction values in different filters as a function 
of A6250- Fortunately, this study avoids the problematic UV region 
where extinction laws ar e characterised by several v ariables, which 
are uncorrected with R jFitzpatrick & Massall2007h . Here, the fo- 
cus is on the IR/optical region where the wavelength dependence 
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of exti ncti on is well defined by R only. Following lSavage & MathisI 
(I1979I) and lHowarthI imA we assumed that = 3. 1, so as to rep- 
resent an average Galactic sight-line. 

The calibration of the broadband e xtinctions to A6250 was 
performed using thelMunarTet "Zl (l2005h library to represent the 
IStraizvs & Kuriliend dlQeil) sequences. Fitting was performed by 
linear regression of a quadratic polynomial to the data, over a range 
of < ^6250 < 10, the functions returned, for a band X, are ex- 
pressed as follows: Ax = axAl^^^ + Z?xA6250 + cx- 

Fig.[6ldemonstrates the effect of extinction on several different 
spectral types in the {r' - i\ r' - Ha) plane. Despite being cosmeti- 
cally similar, there are significant differences in the response to ex- 
tinction of different spectral types. Fig.[7]demonstrates the variation 
of the coefficient for the leading linear term in the fits. For objects 
with types in the critical range A3 to K4 (0.05 - i')< 0.55) 
there is a particularly simple linear relation between this fit coeffi- 
cient and intrinsic {r' -i') colour, as can also be seen in Fig.[71 

The differing responses to rising extinction of different spec- 
tral types is illustrated by the extreme examples of A3V and M4V 
stars: in Fig.[6]these two types limit the plotted {r' - Ha) range of 
the main sequence. As more extinction is applied, the (/ - Ha) 
range spanned shrinks because the A3V star is reddened more 
quickly in this colour. 




(r'-n 

Figure 8. A colour-magnitude diagram for F5-8 stars, selected by 
colour, from IPHAS observations of a 10' x 10' box around {l^b) = 
(101.55, -0.60) . The dwarfs can be seen as the sequence to the left, while 
the giants are seen to the right. Five brighter objects, probably either F 
bright- giants, or highly reddened late-O or early-B stars are visible to the 
extreme right, separated from the giant sequence. 



3 MEAD: AN ALGORITHM TO DETERMINE 
EXTINCTION-DISTANCE RELATIONS 

3.1 Concept 

The basic principle of this work is that we wish to determine the 
distance and intervening extinction to a large sample of objects 
in the IPHAS database. From the observations, we have three ob- 
served parameters for each object, namely the /, i' and Ha mag- 
nitudes, with errors on each. In order to determine the distance to 
an object by photometric parallax, we must determine its absolute 
magnitude and extinction. The simplest way to estimate absolute 
magnitude is to determine the intrinsic colour and luminosity class 
of the object. Thus for each object we have three parameters we 
wish to determine (extinction, luminosity class & intrinsic colour) 
and three observed parameters. 

As described in Section [2l it is possible to ascertain the spec- 
tral class and monochromatic extinction of an object from its po- 
sition on the colour-colour plane. But, there are two main sources 
of degeneracy in this process. The first is caused by the fact that 
reddening tracks for O, B & early A (AO-2) stars intersect those of 
later A and F stars. This follows from the fact that Ha absorption 
peaks in early A stars and (r' — Ha ) is strongly correlated with Ha 
equivalent width JPrew et aDl2005h . As discussed in Section [331 it 
transpires that this is not normally a significant source of degener- 
acy. 

The second more important source of degeneracy is luminos- 
ity class, which has a significant impact on distance determination. 
Therefore, to estimate an object's distance, it is necessary for an es- 
timate of the luminosity class to be determined first. This is, how- 
ever, not a trivial task, as the different luminosity class sequences 
behave very similarly on the IPHAS colour-colour plane. We de- 
cide on the luminosity class of an object as follows. If two ob- 
jects of identical intrinsic colours, but different absolute magni- 
tude, are to have the same apparent magnitude the more luminous 
object must be further away, and is likely to be more heavily red- 
dened. This reddening contrast creates two sequences in a colour- 
magnitude diagram, as can be seen in Fig. [S] where a selection of 
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Figure 9. A schematic depiction of MEAD. 



mainly F stars are clearly split into dwarf and giant sequences. It 
is by identifying these split sequences in the colour-magnitude dia- 
gram that luminosity class is assigned. 
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3.2 Details of implementation 

Fig. [9] outlines an algorithm that realises this principle (named 
MEAD for Mapper of Extinction Against Distance). In the remain- 
der of this section each of the steps in MEAD are described. 

Extinction, intrinsic colour and distance are determined si- 
multaneously for all objects in the intrinsic colour range of inter- 
est given in Figs. |2] and [3] Objects which fall outside this region 
of the colour-colour plane which could be occupied by objects of 
the types we are interested in are simply discarded. The discarded 
objects will include M and late-K type stars, white dwarfs, extra- 
galactic objects and stars exhibiting Ha emission. At the same time 
errors on the derived quantities are also propagated from the pho- 
tometric errors. Although determined at the same time as distance, 
the estimates of extinction and intrinsic colour are almost indepen- 
dent of the luminosity class and can be considered to be determined 
solely from the colour-colour plane: only very small differences 
arise from the slight differences between the luminosity class se- 
quences in the colour-colour plane. 

The appearance of objects on the extinction-distance plane is 
complicated by several factors. Firstly, the size of the errors on the 
estimated values of extinction and distance vary greatly from object 
to object. Also, these errors are correlated with distance and extinc- 
tion, as the apparent magnitude of a source becomes fainter with 
increasing distance or extinction, increased photometric error is to 
be expected. Finally, because of the finite angular resolution of the 
fields, there will inevitably be unresolved substructure within the 
absorbing material in each field analysed, causing an intrinsic and 
irreducible spread of extinctions at any given distance. These com- 
plications are dealt with by assuming that extinction at any given 
distance is modelled by a distribution other than the delta function: 
a maximum likelihood estimator is used to determine the value of 
the parameters that describe the distribution. The choice of family 
of distributions used here is further discussed in Section I441 

IPHAS observations are subject to both bright and faint mag- 
nitude limits. The implication of this is that for an object at a given 
distance, it will only be observed if its extinction also falls between 
two limits, and when these limits are approached, the observations 
will start to become incomplete. This would bias any estimates. 
The nature of the bias is that as the faint cut-off is approached, only 
stars viewed through lower extincti on remain in th e sample, while 
more extinguished objects are lost. iHakkila et al. ([l997.) note the 
disruptive effect that this had on previo us studies ( Fitzgeraldll 19681 : 
Neckel et al. I I198QI :I Arenou et al. 1992). Even at the higher angu- 
lar resolutions treated here this effect continues to be a problem 
because the i nterstellar medium is inhomogeneous down to very 
small scales. iPeshpandd (l2007h demonstrates that variations in the 
structure of the ISM naturally extend down to at least the scale of 
tens of AU. Therefore, it is necessary to continue to attend to the 
magnitude limits in this study, even though significantly finer an- 
gular resolutions are achievable than those in earlier work. 

To deal with this problem, it is necessary to determine which 
objects could be close to the magnitude limit in distance extinction 
space, so that they can be excluded when producing an unbiased 
distance-extinction relationship. To establish which objects are to 
be excluded, first all the objects are binned in distance to produce a 
biased extinction-distance relationship, for each bin the mean and 
intrinsic standard deviation of the extinction is found. Then, for 
each object the maximum (Jmax) and minimum (<imin) distances for 
which observations of that type of object would be complete are es- 
timated. These distances occur where the following two equations 



and faint magnitude limits respectively and o^{d) is the intrinsic 
standard deviation of extinction in the relevant bin: 



A6250('^max)+^a^(<imax) =^6250^ 
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are satisfied, where Ag^f^V) and A^^^^^ 
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If these equations are satisfied at more than one distance, the 
largest possible value of Jjnin and smallest possible value of Jmax 
are taken. In equations [T]and[2]the value of n can be altered. Increas- 
ing it will reduce the chance of the bias affecting the data, at the cost 
of decreasing the number of objects which are binned to produce 
the distance-extinction relationship, whilst decreasing it will act in 
the opposite sense. There is no immediately obvious value to adopt 
for n and Section 14.51 discusses the use of simulated photometry 
to find an appropriate value. At the bright end it is assumed that 
the magnitude range over which observations become incomplete 
is small, relative to the apparent magnitude sampling, so objects 
which are closer than Jmin are simply discarded. In contrast, at the 
faint end where the apparent magnitude sampling is much denser, 
the magnitude range over which observations become incomplete 
is significant. Therefore the contribution of each object to the dis- 
tance and extinction estimates is weighted using a sigmoid function 
as below: 

weight = ^ . (3) 

Where /j is the distance modulus of the object in question, 
A^max is the distance modulus corresponding to Jmax and s is related 
to the distance modulus range over which the sample becomes in- 
complete. A suitable value of s was found empirically to be 3. This 
was achieved by fitting a sigmoid to the faint end of histograms 
of the numbers of objects of a given spectral type against apparent 
magnitude. By applying this weight, objects which lie well beyond 
^^max, where the sample is not complete for all expected extinctions, 
are given approximately zero weight, whilst those well within Jmax 
are given a weight of almost one. 

A Bayesian classification method is employed to estimate the 
luminosity class into which each object falls, given the determined 
distance-extinction relationship, the apparent magnitude of the ob- 
ject, an estimate of the extinction to the object and an estimate of 
its intrinsic colour. The luminosity class sequences are assumed to 
have Gaussian spread in absolute magnitude, with all sequences 
having equal and constant spread. The classification is made by 
evaluating the posterior probability that an object belongs to a lu- 
minosity class, given its extinction, the distance-extinction relation- 
ship and the absolute magnitude of that class. The choice of the 
prior used is discussed in Section |431 

A threshold probability is imposed, such that when the proba- 
bilities of all five classes fall below this, the object is removed from 
consideration in that iteration. This removes objects which do not 
appear to fit into any of the luminosity classes. 

These procedures are put into an iterative structure, as shown 
in Fig.|9] The sense of the iteration can be summarised as follows: 
the distance, extinction and intrinsic colour of each object is es- 
timated; then these determinations are binned to get a distance- 
extinction relation; this is then used to determine the luminosity 
class of the objects; and then the distance, extinction and intrinsic 
colour are estimated again, given the new estimate of luminosity 
class. The process has shown good convergence in all cases exam- 
ined so far. In the range of 5-10 iterations are normally required. 
Convergence is discussed further in Section |4l2l 

MEAD begins with the derivation of the distance-extinction re- 
lation defined by the early A- stars alone, that lie along the lower 
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edge of the main stellar locus in the colour-colour plane. This is 
used to obtain the first estimate of the luminosity classes of all 
objects in the field. It is safe to assume that the great majority of 
objects on this strip are near main sequence, given that early- A gi- 
ants are not only relatively rare but also exhibit slightly weaker 
Ha absorption that lifts them to higher (/ — Ha). This breaks the 
degeneracy between the different luminosity classes and enab les 
distances to be reliably estimated at the outset ( Drew et al. 2008). 

When binning the data by distance MEAD ensures that: the 
depth of the bin is at least 100 pc; there are at least eight objects in 
the bin; the total signal to noise in the bin is at least 130. 

It should be noted that not all available data in the IPHAS 
database are used in this algorithm. Only data flagged as either stel- 
lar or probably stellar in all three bands by the processing pipeline 
are accepted, with the remainder discarded straight away (for more 
details of the data processing see Irwi n & Lewi sl l200ll : iDrew et al.l 
12005: Gonzalez-Solares et al. 2008). As noted already in this sec- 
tion only a restricted stellar intrinsic colour range i s exploited 
(marked in Figs. |2] and [3])- this is due in part to the iHouk et aP 
(Il997h calibration ending at K4V and also to the difficulties asso- 
ciated with modell ing the m olecular bands in late type- stars (Sec- 
tion [J}. Following brew et al. (2008), a maximum bright limit of 
/ = 13.5, = 11.5 and Ha = 12.5 is employed to avoid objects 
that may be saturated in one or more bands. Occasionally when 
any fainter objects are classified as saturated by the pipeline, the 
limit will be moved to the apparent magnitude of that object. 

3.3 Possible contaminants 

MEAD assumes that all the objects examined can be described by 
one of the luminosity class sequences. In practice there are a mi- 
nority of exotic objects that will not fit into them. Generally such 
objects are not a problem if they either: have colours which are 
outside the range MEAD looks at or are significantly less frequent 
than normal objects. The latter argument relies on the statistical na- 
ture of MEAD. Although the properties of individual objects may 
be incorrectly estimated, the intention is that MEAD will derive the 
correct extinction-distance relationship given a large enough sam- 
ple. 

As previously mentioned, O and B stars are contaminant ob- 
jects that may be mistaken for A and F stars. Fortunately, O and 
B stars are relatively rare, and present only in young populations. 
In populations older than ~100 Myrs they are essentially absent. 
So, in large samples, with a mix of ages, we can comfortably as- 
sume field A and F stars dominate. The circumstance in which this 
assumption may come undone is in regions of ongoing star for- 
mation, especially if the stars that will become A and later-type 
main sequence stars are still in the pre-main- sequence phase. If 
such objects remain in the minority for their locale, MEAD should 
recognise their raised luminosity and classify accordingly. If they 
become the dominant population, then MEAD may misinterpret 
them. In Section [53l MEAD is tested against IPHAS observations of 
the Cygnus OB 2 association, where such extreme conditions could 
exist. 

The subdwarf sequence occupies a very similar position on 
the IPHAS colour-colour plane as its dwarf equivalent and as such 
there is the possibility that subdwarfs may be erroneously identified 
as dwarfs and used in MEAD. However, in a volume limited sample 
it is expec ted that roughly one in every 10-^ objects would be a 
subdwarf dReid etalJl200lh . Given the faint absolute magnitudes 
of subdwarfs relative to the dwarfs, this proportion becomes even 
lower in the magnitude limited IPHAS catalogue. Furthermore, the 



imposition of a threshold probability in the Bayesian luminosity 
classifier will prompt MEAD to discard a large proportion of them 
as they are offset below the main sequence in colour-magnitude 
space. 

Hydrogen atmosphere white dwarfs largely lie outside the 
main stellar locus, leaving only a minority inside the region of the 
colour-colour plane searched by MEAD: these objects will be the 
very coolest and the very hottest white dwarfs. White dwarfs ex- 
hibit absolute magnitudes approximately in the range 10 < / < 15. 
Given this and the IPHAS magnitude limits, the database cannot 
contain even unreddened white dwarfs more distant than 1 kpc. 
It is estimated that there will be 7.6 white dwarfs per square de- 
gree in the IPHAS sample, representing ~ 0.01% of all sources. 
This estimate is based on a detailed mod el of the Galactic popu- 
lation of white dwarfs dNapiwotzkil 2008h calibrate d us ing the ob- 
served white dwarf samples o f'Holberg et al.72008') and P auli et al.l 
12006), reddening is included in an approximate way. Many white 
dwarfs will be ignored by MEAD due to their colours lying outside 
the main stellar locus. Of the handful remaining within more will 
be discarded by the Bayesian luminosity classifier as belonging to 
a faint sequence. 

Extra-galactic sources will also contaminate the catalogue. 
Resolved galaxies will be flagged by the CASU pipeline as be- 
ing extended objects and will be ignored by MEAD. In addition, 
even some small galaxy images will overlap with the very numer- 
ous stellar images and hence be classified as extended, also remov- 
ing them from consideration. Nevertheless, some galaxy images 
will be sufficiently compact to be classified as stellar. Genuinely 
compact (i.e. physically small) galaxies are rare, contributing about 
3% of all l ow redshift field galaxies down to 5 ^ 20 (/ ~ 19) ac- 
cording to Drinkwate r et al.l (Il999|). who used the all-ob ject (high 
latitude) spectroscopic survey of iDrinkwater et al ] (l2000n . Extrap- 
olating this another magnitude deeper (i.e. to / = 20), we would 
still only expect some 40 objects per square degree, which is com- 
pletely negligible compared to (i.e. 0.02% of) the 200000 stars per 
square degree in low latitude fields. Of course, at faint limits all 
galaxies tend to look less and less resolved, resulting in the well 
known merging of the st ellar and galaxy se quences in star/galaxy 
separation diagrams (e.g. [Driver et al ]|2003V However, at / = 20 
there are around 2000 galaxies per square in total, a pproximately 
two th irds of which are in the final magnitude bin (e.g. [Driver et al.! 
Il994l) . so even if we assume that all these could look stellar, this 
would only represent a 0.7% contamination of the "stars". In re- 
gions of high extinction, the counts would, of course, be reduced at 
given /. For an asymptotic Ar = 2, for instance, the galaxy counts 
would be reduced by a factor of about 8, though this would be com- 
pensated somewhat by the galaxies having smaller angular sizes 
as the low surface brightness outer parts are too faint to detect 
(Phillipps. Ellis & Strond [T981 ). This would reduce the contami- 
nation to ~ 200 per square degree or ~ 0.1%. 

Quasars and other distant AGN will also show unresolved 
images, but again numbers ar e small . From the IDrinkwater et aD 
(2000) all-object survey, Mev er et al.l (1200 ih find 35 QSOs (of all 
colours) per square de gree to ^ 20.2 (/ ^ 19.5). Given the steep 
quasar number counts JCroom et al.ll2004l), th i s wou ld be ~ 70 per 
square degree at / = 20. IVanden Berk et al.l (l2005h find a some- 
what lower number density in the band and many QSOs may 
have colours too blue to be included in our sample here, so we can 
estimate a contamination rate of order 0.02% (and less, of course 
in highly absorbed regions). 

Finally we note there are several classes of emission line star, 
including symbiotic stars, cataclysmic variables and planetary neb- 
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Figure 10. The growth of random photometric errors in real data (red), from Figure 12. The determined points on the extinction-distance relationship 

IPHAS field 231 lo and simulated data (black), given similar seeing. from twenty five simulations (black) and the underlying extinction, input to 

the simulations (red). 



ulae, which lie above the main stellar lo cus on an IPHAS c olour- 
colour plot and so are ignored by MEAD jCorradi et alJl2008h . 



4 TESTING MEAD AGAINST SIMULATED 
PHOTOMETRY 

The reliability and precision of MEAD needs to be appraised, by 
testing its performance on sightlines where the distance-extinction 
relation is exactly known. As there are no such ideal sightlines, it 
is necessary to synthesise model sightlines resembling those en- 
countered in the Galactic disc, to which photometric errors can 
be applied that are consistent with IPHAS data. We do this using 
Monte-Carlo sampling techniques. These simulations also allow us 
to directly investigate the impact of factors such as metallicity on 
the performance of MEAD, in a way that would otherwise not be 
possible. 

The simulations were created largely following the Besangon 
galaxy mode l. The IMF of [K roupa (2001) and the galactic extinc- 
tion model of lAmores & Lepine. C2005) were assumed. The simula- 
tions were populated with objects from a copy of t he Teramo library 
of isochrones and stellar evolutionary tracks (Pietrinferni et al. 
12004 ). These were translated into the IPHA S filter system usin g 
the finely-sampled grid of model SEDs due to lMunari et al.l (l2005h . 
Photometric errors were imposed, adopting random errors based 
on the CCD equation (see Fig.[TO]) along with field wide systematic 
errors of order 1.5%, arising from the calculation of the aperture 
correction and photometric zero-points. 



4.1 Comparisons between simulations and real data 

Although it would be possible to create an end to end simulation, 
whereby synthetic images are created and then pipelined in the 
same manner as IPHAS observations, this would be so excessively 
time consuming as to restrict the analysis of mead's performance. 
Instead, the simpler method used here begins with the application 
of MEAD to the synthesised magnitudes. This allows the production 
of many synthetic lines of sight, which in turn enables MEAD to be 
robustly analysed. 

The colour-colour and colour-magnitude plots produced by 
the simulations bear strong similarity to the observed plots, as 



shown in Fig. \TT\ though there are a few key differences. Signifi- 
cantly, the clear group of M-dwarfs are not reproduced in the simu- 
lations, due to the lower mass limit of the Teramo library. This is not 
a problem in this context, as MEAD does not use objects which oc- 
cupy the region of the colour-colour plot occupied by M-type stars. 
The same argument applies to the exotic object types discussed ear- 
lier (see Section [331) . These are not included in the simulations, as 
they do not significantly influence MEAD. 

The difference between the two plots portrays the divergence 
that is to be expected between Galaxy models and the details of 
the real Galaxy. The Amores & Lepine ( 2005) and real extinction- 
distance relations are in clear disagreement, giving rise to the 
different angles made by the main stellar loci with respect to 
the axes in the colour-magnitude diagrams of Fig. \TT\ The real 
colour-colour diagram in Fig. \TT\ shows two principal groupings, 
one with low reddening, as seen in the group that runs through 
(/ — — Ha) ~ (0.6,0.3), while a second group with higher 
reddening occurs around ^ (0.9,0.35). These two groupings occur 
because the stars in this direction are not distributed smoothly with 
extinction, but rather are either local unreddened stars or are more 
heavily reddened and therefore probably more distant. The discon- 
tinuity may arise either as a result of the stars clumping together 
into two groups at different distances, or as a result of a sudden in- 
crease in extinction. The simulations, on the other hand, do assume 
that stars are distributed smoothly and extinction grows steadily, so 
there is no grouping visible. 

4.2 The basic performance of MEAD 

With the photometric errors, metallicity variation and binary frac- 
tion set to zero, MEAD is very successful at retrieving the input 
distance-extinction relation, as shown in Fig. [TJl The little spread 
on extinction that does exist is a result of describing the continuous 
stellar evolutionary tracks in terms of 5 discrete luminosity classes 
(see Fig. [3 . In this manner we can see that the errors induced 
by this description are small (5A^/ ~ 0.05) and so this is a valid 
method. 

To be credible, the solution must converge. Although MEAD 
will exit when complete convergence has been reached it is unrea- 
sonable to wait for every last small change to be made, especially 
if the last corrections to a few stars become oscillatory. Fig. [13] 
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Figure 11. A comparison between observed (left) and simulated (right), colour-colour (top) and colour-magnitude (bottom) diagrams , given similar conditions, 
in the direction {l^b) = (170,0). The observed data are the observations of IPHAS field 231 lo and include only objects classified as stellar or probably stellar. 
A magnitude cut at r' = 20 has been applied to both the observed and real data. The differences between the two sets of data are due to: differences between 
the real and model distance-extinction relationship; differences between the distributions of objects with respect to distance; and some types of objects not 
being included in the simulations. 



demonstrates the convergence rate for one execution of MEAD, 
whilst Fig. [T4l characterises the convergence rate for a large num- 
ber of visualisations. The default maximum number of iterations 
for MEAD is set to 15. 



4.3 Calibration of priors for the determination of luminosity 
class 

In Section [3] it was noted that it was possible to include a non- 
uniform prior probability when determining the luminosity class of 
a given object. Using the simulated photometry the probability dis- 
tribution for the prior which retrieves the most accurate result was 
determined from a representative sample of sightlines. The accu- 
racy is defined in terms of minimising the statistic: 



where Ag^^ is the extinction estimated by MEAD and A^^? the extinc- 
tion from the model at the same distance. In the denominator, 
is the intrinsic standard deviation of extinction in the distance bin 
(i.e. due to the porosity of the ISM), while 5A/ is the error in Aest 
as a consequence of the photometric errors. 



Given that simulating sufficient sightlines is relatively time- 
consuming, it was not possible to search the entire parameter 
space, whilst the uncertain structure of the parameter space made 
'downhill-only' root finding algorithms inadvisable. Therefore, a 
simulated annealing algorithm was used to minimise the statis- 
tic above. It was found that the prior should take the following 
form, where P{X) is the probability assigned to luminosity class 
X: Piy) = 0.34; P(/y) = 0.31;P(///) = 0.28; P(//) = 0.06; P(/) = 
0.01;. This prior mirrors the number of objects in each class in 
the simulated photometry. As the IPHAS bright magnitude limit 
is at / ~ 13.5, neither bright-giants (M^/ < —2) nor super-giants 
(Mr' < —6.5) are likely to be present in large numbers. 

We note that this set of priors has been derived from a weight- 
ing of sightlines that, at the present time, emphasises the moder- 
ately reddened directions in the Galactic Plane encountered at all 
latitudes outside the Solar Circle, and at latitudes |^| > 1 within it. 
This set was specifically chosen to correspond to those sightlines 
investigated in this work. For these sightlines dwarfs are strongly 
dominant and so the results show little sensitivity to the choice of 
priors. The priors for the highest reddening sightlines passing into 
the inner Galaxy are under investigation. It may prove necessary 
to attribute somewhat increased probabilities to brighter luminosity 
classes. 
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Figure 13. Demonstrating the convergence of MEAD by showing how the 
number of objects which have their estimated luminosity class changed 
varies with iteration. From iteration 10 onwards MEAD is oscillating be- 
tween two solutions, in doing so the luminosity class estimates of four ob- 
jects are changed in each iteration. 



Figure 15. The difference between the values of mean monochromatic ex- 
tinction at 6250A at a given distance determined by MEAD (Aout) and the 
value in the simulation (Ain). The simulation has had intrinsic scatter in ex- 
tinction imposed in a manner consistent with Kolmogorov turbulence in the 
ISM. 
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Figure 14. 1^ is the number of iterations required for the number of objects 
which have their estimate of luminosity class altered to drop by a factor of 
e, as determined by linear regression. Ig is a measure of the gradient of plots 
like Fig. [13] where lower values of Ig correspond to steeper gradients and 
so quicker convergence. 



Gaussian, then the estimate of the mean value of the extinction in 
each bin (A) is given by the following equation: 

y M 

where A/ is the estimated extinction for a single object, and other 
terms are as described in Section |431 

This estimate is essentially a weighted mean of the individual 
estimates of extinctions for each object in the bin. It is found that 
the intrinsic scatter of extinction in each bin (a^) can be satisfacto- 
rily estimated by the weighted standard deviation of the extinctions 
of the brighter objects in the bin (those with dA^25Q ^ O-l)- The 
weighting for this estimate is the weight given by equation [5] over 

In reality, the intrinsic scatter of extinction in each distance 
bin is a result of the turbulent processes at work in the ISM. If the 
extinction of each simulated object is scattered away from the mean 
value for the distance as a result of Kolmogorov turbulence, MEAD 
still successfully retrieves the input extinction, despite the in-built 
assumption that extinction at a given distance is scattered normally 
(Fig.CSJ. 



4.4 Coping with unresolved substructure & photometric 
errors 

In MEAD, once estimates of distance and extinction are produced 
for each object, the data are binned, with respect to distance, to pro- 
duce a distance-extinction relationship. As discussed in Section [3] 
the extinction will not be uniform for a finite sized box at any given 
distance. Instead it will vary with sky angle, just as it will vary 
with depth, due to unresolved substructure prompting some intrin- 
sic variation in extinction. Hence, in this study, the extinction in 
any given bin is characterised as having some intrinsic scatter. 

If the photometric errors are propagated through to give the 
error on the determined value of the extinction, it is possible to 
determine Maximum Likelihood Estimators (MLEs) for the mean 
and intrinsic standard deviation of the extinction, if the character of 
the variation in extinction is known. If the variation of extinction is 



4.5 Handling of magnitude limits 

In Section [T2l we discuss the difficulties associated with the magni- 
tude limits and how these have affected previous extinction maps. 
The treatment for this problem involves ignoring data which are 
at a distance where they could be affected by incompleteness. The 
parameter n in equations [T] and |2] determines which data are ig- 
nored, but there is no obvious value it should take. To establish the 
ideal value of «, MEAD was tested against variation of this param- 
eter. Setting the value of n too low allows the inclusion of objects 
which may bias the results, due to the effect of incompleteness be- 
coming increasingly relevant. It was found that for ;i < 2 the re- 
sults of MEAD would be unsatisfactory, an example of this effect is 
shown in Fig. [16] Raising n acts to reduce the number of objects 
which can be used to determine the distance-extinction relation- 
ship, thus decreasing the range and distance resolution of the out- 
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Figure 16. Setting the value of n too low in equations [T] and |2] thus 
compromising the results of MEAD. Here the upper panel shows n = 0.5 
and the lower n — 2.0. Both panels are for 50 iterations in the direction 
{l^b) = (180,0), the actual (input) distance-extinction relation is the red 
line and the black points are the results from MEAD. 




Figure 17. The variation in the total number of data points (black dots, 
left axis) and proportion of bad data points (red crosses with error bars, 
right axis) on the derived distance-extinction relation as a function of n for 
the directions {l,b) = (80,0), (130,0), (180,0). Where good bins are those 
which satisfy the condition Agut ~ 

< + 6^ and bad bins are those 

which do not. 



come. Although many of the objects filtered out should be removed 
for reasons of incompleteness, others need not be and are lost un- 
deservedly. A compromise is necessary. We set = 2 so as to use 
the maximum possible number of objects whilst avoiding the mag- 
nitude limit effect. Fig. [TT] illustrates that setting n = 2 obtains the 
highest number of accurate distance bins on the derived distance- 
extinction relationship whilst maintaining a low proportion of bad 
bins. 



4.6 Achievable angular resolution 

Previous observational 3D extinction maps of the Milky 
Way (Fitzgerald 1968; Lucke 197^ iNeckel et all Il98d : 
IPandev & Mahral 119871 : iJoshil 120051 . e.g.) have been carried 
out at angular resolut ions of the order of a few degrees. The recent 
extinction map due to [Marshall et has a resolution of 30 . 

With the high angular density of sources in the IPHAS database 
(~ IQr^ sq.degr^) it is clearly possible to map extinction at a much 
finer angular resolution. The resolutions possible are dependent on 
the source density, but even in the relatively sparsely populated 
galactic anti-centre resolutions of 10' are achievable. 



4.7 Bias in distance estimation 

The estimate of the distance modulus to each object has a symmet- 
ric error distribution around its mean. However, the expected distri- 
bution of stars is not symmetric as the volume of space in which a 
star lies and the density of stars changes. Therefore, the expectation 
of the true distance to the star does not equal the original estimate 
and so that original estimate is biase d. An analogous effe ct leads 
to the well known Lutz-Kelker bias (iLutz & Kelkei|[l97 3S. when 
using trigonometric parallaxes to esti mate t he absolute magnitudes 
of objects. Following Lutz & Kelkerl (Il973h it is possible to derive 
a correction for this bias, if the density of stars is assumed to be 
uniform. If this correction is applied to the distance bins of ob- 
jects from the simulated photometry, where the stellar density is 
assumed to be uniform, the median correction is found to be 0.1%, 
whilst 95% of bins have corrections less than 1%. The improvement 
on the agreement between the input and derived extinction-distance 
relationships is thus negligible. 

Calculating a correction for a realistic density distribution is 
considerably more difficult. However, the median correction will 
often be smaller than in the case of uniform stellar density, because 
of the drop off in stellar density associated with most sightlines. 



4.8 Sensitivity to disc metallicity gradients 

The Besan9on model assumes a smooth gradient in the mean 
metallicity with resp ect to ga la ctocen tric ra dius {Rg). This 
view is shared by iFriel et all (l2002h and 'Bragad ia & Tosil 
(2006) among st others. However, i t has been ar gued (e.g. 
Twarog. Ashman~& Anthonv-Twarogj 1 19971: [ Corder & Twaro j 



200 ll : lYong. Carnev & Teixera de Almeidal l2005l) that metallicity 
in fact varies as a step functi o n of R g, with the step oc c urring 
at Rg ^ 10 kpc. ICarraro et all tOOf i and "Bragaglia et aL| (l2008h 
argue that any gradient flattens to a constant metallicity in the 
outer disc. Given this spread in opinion it was necessary to test the 
performance of MEAD against a range of different schemes. 

It was previously shown in Fig. |4] that metallicity has almost 
no effect on the colours of A0-K4 stars. The absolute magnitudes. 
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on the other hand, are metalUcity dependent, but these do not be- 
come significant until metalUcities are markedly sub-solar (Fig. [5]). 
Consistent with this we find that introducing metallicity variations 
induces little change in the output from MEAD, even though MEAD 
assumes solar metallicity. This is true of all the above mentioned 
schemes for the disc metallicity gradient. 



4.9 Correction for binarity 

Unresolved binaries and higher multiples present the possibility 
of the systematic over-estimation of point- source absolute mag- 
nitudes. For the objects selected by ME AD (A0-K4), the effec t of 
binarity on the observed colours is small (iHurlev&Toutll 19981) . 

In the extreme case of binaries of unit mass ratio the distance 
to an object would be underestimated by a factor of \/2 compared 
to the single- star case. However, on making plausible allowances 
for the binary mass ratio distribution and for the binary fraction, 
it is likely that, in the mean, distances will only be underesti- 
mated by a factor of a few percent. For exa mple, the expected un- 
derestimate of distance will be ~ 5% if the iKroupal 42001) IMF, 
iDuquennov & Mavo3 (Il99lh binary fraction, a constant probabil- 
ity function for binary mass ratio and a mass-luminosity relation of 
L oc M^-^ are all assumed. 

To arrive at a correction factor that may be applied routinely 
within MEAD, the effect of binarity has been analysed through the 
use of simulated photometry. Binaries were inserted in the model 
using a binary fraction of 57%, following IPuquenn ov & Mavoj 



19911) which is derived from observations of G-type stars. Lada 



20061) demonstrates that the binary fraction drops significantly 
for later type stars, but these are mostly those which are not 
simulated (i.e. M-type stars and later). The probability function 
for the binary mass ratio was assumed to be constant, whilst 
higher order unresolved mul tiples are assumed to be relatively rare 
touquennov & Mavo3ll99lh and so are not included in the simula- 
tion. 

The outcome of this exercise is that there is systematic un- 
derestimation of distances, which in turn causes extinction to be 
overestimated in a given distance bin. The appropriate correction 
for this is made by multiplying all the derived distances by 1.06. 
This figure was derived through minimisation of the statistic in 
equation m over 100 visualisations of 6 different sightlines. This 
modification should function successfully as most bins contain rel- 
atively large numbers of objects and the correction itself is small. 
When this correction is adopted MEAD successfully reproduces the 
input distance-extinction relationship. 



4.10 Impact of reddening law variations 

MEAD and the simulation s performed up unt il this point assume the 
R = 3.1 reddening law of lFitzpatrickl (|2004|) . To investigate the ef- 
fect of altering R significantly, simulated photometry was created 
with R = 2.6 and 3.6 reddening laws. Then MEAD was run on the 
results of these simulations, while continuing to apply the = 3.1 
law. Fig.fTSlshows that the impact on MEAD's mapping of sightlines 
better represented by non-standard reddening laws is to induce a 
systematic error that is comparable - in these examples - with the 
typical random error. Had the reddening been expressed in terms 
of the monochromatic extinction, Fig. [18] would look much 

the same. However, a conversion into Ay, would be even more de- 
pendent on the value of R and therefore the outcome of a mapping 
expressed in terms of Ay would be subject to even greater system- 
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Figure 18. The effect of altering the value of R on the difference between 
the values of mean broadband extinction in the r' band at a given distance 
determined by MEAD {A{r')out) and the value in the simulation {A{r')in). 
The results for = 3.6 are represented by unfilled black bars bounded by 
solid lines, = 2.6 by hatched red bars bounded by dashed lines. The sim- 
ulations are for the direction (l^b) = (180,0) and 50 visualisations were 
performed in each case. 



atic error. The sense of the bias is that the extinction at a given 
distance for sightlines of anomalously high R is underestimated. 

Given that IPHAS observations provide magnitudes in just 
three bands, it is not feasible to constrain R in addition to the in- 
trinsic colour, extinction and luminosity class (and so distance). 
Instead it is preferable to rely on the relative dominance of 
the = 3.1 reddening law across the Gala ctic Plane. The sarn e 
assu mption has been made previous l y by iNeckel et al.l (1980'); 
Schleg el et all (Il998h : [Marshall et aP (l2006h and many others. If 
the IPHAS database were to be cross-matched with the results of 
the United Kingdom Infrared Deep Sky Survey Galactic Plane Sur- 
vey (UKIDSS GPS: lLucas et al.ll2008l) and the Ultraviolet Excess 
Survey of the northern Galactic Plane (UVEX; Groot et al. in prepa- 
ration) the extra measurements would facilitate the accurate esti- 
mation of R, whilst still maintaining the accuracy of the other esti- 
mates. 



4.11 The impact of very young populations 

O, B and very early A stars occupy the same area of the IPHAS 
colour-colour plane as less-extinguished later A and F stars. How- 
ever, they also respond differently to extinction and have different 
absolute magnitudes and so introduce a degeneracy into the prob- 
lem. In the inter- arm field such objects are comparatively rare, with 
their short lifespans and the IMF suppressing their frequency. How- 
ever, in the spiral arms such objects are more frequent, due to the 
on-going star formation induced by the passage of the spiral density 
wave. The Besangon model makes no allowance for spiral-arm in- 
duced star formation, since a constant star formation rate is adopted 
for all locations. Therefore simulations based on the Galactic model 
cannot provide any insight into the impact of these localised young 
populations. We have carried out some tests of such concentrations 
by inserting them into model sightlines in order to mimic the super- 
position of a spiral arm. The results show no additional systematic 
error in the derived extinction curves that would imply a breakdown 
of the method. 

As with early- type stars, pre-main sequence (PMS) objects in- 
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troduce degeneracy, due to their different colours, absolute mag- 
nitudes and response to extinction as compared to the sequences 
assumed in the calibration of MEAD. Also, as with early-type stars, 
PMS objects would normally be found concentrated in star forming 
regions. So as to investigate the effect of including PMS objects the 
extra spiral arm population was again included in the model. The 
inserted PMS objects were then tracked through MEAD. In reality, 
many PMS objects exhibit Ha emission and so lie outside the main 
stellar locus on colour-colour plots and are thus ignored by MEAD. 
However, those PMS objects exhibiting weaker Ha emission would 
remain inside the main stellar locus, where they may be misiden- 
tified by MEAD. As with O and B stars, the relative scarcity of 
emission-line PMS stars compared to normal stars should limit this 
problem. As a direct test of this and other expectations regarding 
young populations, we present and discuss the example of sight- 
lines through Cyg 0B2, a spectacular northern OB association, in 
section 1531 



4.12 The impact of photometric zero-point errors 

At the time of writing the IPHAS database lacks a global calibra- 
tion, with observations instead being calibrated on a run by run 
basis, with the Ha calibration tied to the r' calibration. As such 
it is possible that there will be offsets between the true and mea- 
sured magnitudes in one or more bands. To investigate the effects 
of such difficulties, offsets were manually introduced into the sim- 
ulated photometry. 

The effect of adding the same magnitude offset to all three 
bands is simple to understand, as this will simply lead to inaccurate 
determinations of the distance to each object. For example an offset 
of 0. 10 mags is well within the range of what might be encountered 
within the IPHAS database presently. This offset leads to a 5% error 
in determined distances. 

Somewhat less straightforward is the effect of offsetting one 
band only, as this will affect the intrinsic colours and extinctions de- 
termined for each object, and so the object distance as well. Fig. [19] 
shows the effect of applying a positive and negative offset to Ha. 
There are two consequences of doing this. The first is that the ac- 
curacy of MEAD is reduced, as all objects are now systematically 
misidentified. This effect can be seen clearly by comparing the two 
panels in Fig. \\9\ Second, MEAD returns fewer points, for the case 
of a negative offset to Ha. This can be very serious as it may re- 
duce the number of returned points significantly. In the case shown 
90% of the objects available are lost: this is because the main stel- 
lar locus has shifted largely outside the region of the colour-colour 
diagram searched by MEAD. 

Until a global photometric calibration is performed, tying the 
Ha calibration to / vastly reduces the frequency and size of offsets 
in the {/ — Ha). In the future the calibration in the band will also 
be tied to the / band calibration, improving the suitability of the 
data for use in MEAD. 



5 EXAMPLE APPLICATIONS OF MEAD 

5.1 A sightline in Aquila, through (l,b = 32.0, 2.0) 

For this direction, the results from this work (Fig. (20]) show two 
significant increases of extinction, one within the first few hundred 
parsecs, which would be associate d with the Aquila Rift. The Rift is 
at a distance of ~ 200 pc ( Dame et al.ll987), with the extinctio n ris- 
ing to Ay ^ 3 (Ar' ^ 2.5) dStraizvs. Cernis & Bartas"iat3l2003h . We 
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Figure 19. The top panel shows the retrieved extinction-distance relation- 
ship, given an offset of 0.05 (black dots) or —0.05 (green crosses) to the 
Ha magnitudes of the objects. The bottom panel shows the results with no 
offset (black dots). On both panels the input relationship is shown as the red 
line. The simulations are for the direction (/,/?) = (180,0). 




T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r- 



4 6 
Distance / Icpc 



10 



Figure 20. Distance-extinction relationship for field 4199 in the direction 
{l,b) = (32.0,2.0), using a 10' x 10' box. The relationship derived in this 
work is in black, with the error bars indicated the e rrors on the derived 
mean value of extinction in each bi n. The result of iNeckel et al.l (ll98Qh 
is in green, the ISchlegel et al.l (Il998|) asymptotic value is shown with the 
blue dashed li ne, the [Primmel & SpergeJ bOOlh model in cyan and the 
Marshall et al] ( l2006h results are shown with red crosses. Estimates of ex- 
tinction and distance for 1197 objects have been binned up to produce the 
distance-extinction relationship. 
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Figure 21. The retrieved extinction-distance relationship for (/,/?) = 
(177.6,3.1) (field 2849), for a 10' x 10' box. The distance-extinction re- 
lationship derived in this work is in black, with the error bars indicated the 
errors on the derived mean value of extinction in each bin. The result of 
iNeckel et"S1 17980) is in green, the lSchlegel et al.l Il998 h asymptotic value 
is shown with the blue dashed line and t he distance and extinction of NGC 
2099 as found by 'H artman et al.l 1200 8^ is shown in red. Estimates of ex- 
tinction and distance for 752 objects have been binned up to produce the 
distance-extinction relationship. 



find that the extinction through the Rift is somewhat higher, reach- 
ing Ay ~ 4.7 (or Ayi ~ 4) at this position. The second extinction 
rise lies at ~ 2 kpc and is likely to be associated with the Sagittarius 
Arm. The extinctio n then asymptotes be yond it at a value roughly 
consistent with the Schlegel et al. fl998) limit for this direction. 

Previous extinction-distance relations for this sightline are 
also shown in Fig. |20l At distances greater than 2 kpc, MEAD pro - 
duces a relation that is quite simil ar to that of^Marshall et al.l (l2006h . 
It is clear, however, that the Mars hall et al.l (l2006h modelling does 
not take into account the more local Aquila Rift. It is this differ- 
ence, in particular, that lies behind the great discrepancy between 
their first data point at ~2 kpc and our relation. At all distances we 
find the extinction to be markedly higher than given by the relation 
due tolprimmel & Spergel (2001). These distinctions begin to il- 
lustrate the value of an approach that does not depend on Galactic 
models. 



5.2 The sightline through the open duster NGC 2099 

177.6,3.1) 

The Galactic cluster, N GC 2099 (M37) lies near the Anticentre 
at {l,b) = (177.6,3.1)). iHartman et aP (l2008h measure this clus- 
ter to be at a distance of 1.49 ±0.12 kpc with a reddening of 
E(B -V) = 0.23 ± 0.04. This sightline falls within IPHAS field 
2849. The distance-extinction relationship derived in this work and 
in the literature, are shown in Figures [21] and |22] The difference 
between the two plots is that in the first the errors on the esti- 
mated value of extinction for each bin are marked (these arise from 
systematic and random errors involved in obtaining the photome- 
try), while in the second, the error bars reflect the intrinsic scatter 
that is a consequence of the existence of unresolved substructure in 
the interstellar medium. Both of these must be known to estimate 
the distance of an object based on its extinction and the distance- 
extinction relationshi p. Using the distance -extinction relationship 
shown above and the lHartman et al.l (l2008h estimate of reddening, 
the distance to NGC 2099 is estimated to be 2.0 ±0.9 kpc, where 
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Figure 22. As with Fig. [21] only here the error bars indicate the estimated 
intrinsic variation of extinction in the bin. 
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Figure 23. Distance-extinction relationships for the IPHAS field 5985, in 
Cyg OB 2, which has been divided into four quadrants. The relationship 
for the south-east quadrant is indicated by black dots, the south-west by 
red crosses, the north east by green open diamonds and the north west by 
blue triangles. Also shown are the mean distance extinction curves for fields 
6010 (black solid line) and 6134 (red dashed line). 



the large error is a result of the extinction-distance relationship flat- 
tening off between ~ 2 kpc and ^ 3 kpc. Based on the distance to 
W30H of 1.95 ± 0.04 kpc (IXu et al.ll2006 ). this void in extinction 
wou ld appear to lie directly b ehind the Perseus spiral arm. 

lNilakshi&Sagaij(l2002h measure a reddening of E{B — V) = 
0.30 ±0.04 and a distance of 1.36 ±0.1 kpc to NGC 2099, this 
estimate is a significantly poor er fit to derived dist ance-extinction 
relation than the estimate of .Hartman et al The cause of 

the differen ce between these two results appears to be the adopted 
metallicity: Hartman et al. (2008') adopt a nearly solar value, while 
Nilakshi & Sagar ( 2002) prefer a reduced metallicity of Z = 0.008. 
If all objects we re assumed to have fa inter absolute magnitudes, in 
line with the .Nilakshi & Sa^ai (I2OO2I) metallicity, the distance and 
extinction determined bv lHartman et al.l (^008) remains a better fit 
to the produced distance-extinction relation. 
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Figure 24. The ratio of the numbers of near Main Sequence early- A stars 
to F dwarfs as a function of distance for the field 5985. 



5.3 Cygnus OB2 

A potentially difficult environment for MEAD to operate in is one 
with substantial on-going star formation (Section R.lll ). In such re- 
gions there will be relatively large numbers of pre-main sequence 
objects which are not explicitly identified in the algorithm and 
therefore may be more subject to misclassification, and in addi- 
tion any O and B stars present will be mistaken for A/F stars. One 
such region is the Cygnus 0B2 association. Cy g OB 2 has been ex- 
tensively studied in the past, and in particular ' Drew et aP (l2008h 
utilised IPHAS photometry to study the positions of early A-stars 
around the association, concluding that it lies at a distance between 
1.45 kpc and 1.75 kpc, with the range in distances a result of un- 
certainty in the age of the association. 

IPHAS fields 5985 and 6010 look towards the approximate 
centre of the Cygnus 0B2 association, having field centres at 
= (80.3,1.2) and {l,b) = (80.2,0.7) respectively. Figure|23] 
show the derived distance-extinction relationships for these fields: 
in both cases a substantial increase in extinction is apparent at 
~ 1.5 kpc. This rise in extinction would therefore appear to be as- 
sociated with Cygnus 0B2 association. In field 5985 the extinc- 
tion flattens out and approaches its asymptotic value. This flatten- 
ing is not seen in the results for field 6010 - most likely due to the 
extinction continuing to rise significantly beyond Cyg OB 2, plac- 
ing most objec ts outside the magnitude limits of the observations. 
ISchlegel et all ([1998) find significantly higher values of asymptotic 
extinction in field 6010 than in 5985, consistent with this (see also 
iDrew et al.ll2008l) . 

The range over which the extinction sharply rises in field 5985 
is ~ 400 pc. It is not likely that the extinction associated with the 
Cygnus 0B2 association itself is this deep, but rather that the cu- 
mulative errors are inducing a smoothi ng of the distance -extinction 
relationship (a point commented on bv lJuric et alj|2008 h. Using the 
Monte-Carlo Method for simulating photometry, it is possible to 
produce a result similar to the ones displayed in both fields with 
the extra extinction confined to a depth of < 100 pc at a distance of 
(1.5 ±0.2) kpc in the model. 

A second method for determining the distance to the Cygnus 
0B2 association is to determine the ratio of the number of younger 
stars to a sample with longer lifetimes along the line of sight. Early 
A (AO-5) near Main Sequence stars were used as the younger pop- 
ulation due to their relative ease of selection and the low probabil- 
ity of contamination by non-MS stars. Their number was compared 



with that of late F-type dwarfs (F5-9), a population with a main se- 
quence lifetimes roughly 10 times longer. Ideally a sample of red- 
der dwarfs or giants would be preferred in this role, but later-type 
dwarfs are too faint to be observed over a large enough distance 
range, due to the rapidly-rising extinction, while giants are not suf- 
ficiently numerous in this region. Fig. l24l demonstrates the depen- 
dence of this ratio on distance within field 5985, showing a signif- 
icant peak at 1.55 kpc. This distance is consistent with both that 
determined in Drew et al. (2008) and the rise in extinction shown 
in Fig. [23] The peak shown in Fig.[24lis not as broad as the region 
of high differential extinction in Fig. (23] this is a consequence of 
the near Main Sequence early- A stars being relatively bright and so 
exhibiting smaller photometric and therefore distance, errors. 

As previously mentioned, a significant concern when 
analysing fields in Cyg OB 2 is contamination by PMS objects. 
PMS objects will exhibit brighter absolute magnitudes than MS 
objects of the same intrinsic (/ — /') colour and thus their distance 
would be underestimated if they were classified as dwarfs. There- 
fore it might be expected that distance bins closer than Cyg OB 2 
might incorrectly contain objects with higher extinctions, thus 
causing the estimate of mean extinction in these bins to be too 
high. However, many PMS objects will actually be classified in 
luminosity classes IV or III, especially if they have blue intrinsic 
colours where the luminosity classes are tightly spaced (Fig. [2]). 
The ability to classify these objects in classes IV or III will result 
in a more accurate estimation of distance. But it must be acknowl- 
edged that there is room for some falsification through application 
of absolute magnitude scales derived for evolved stars. A second, 
potentially more serious, problem is a result of the behaviour of 
PMS objects on the colour-colour diagram, particularly in relation 
to the Ha band. Although substantial Ha line emission will move 
an object outside the main stellar locus, less pronounced emission 
could result in an object still lying within the main stellar locus and 
MEAD estimating an intrinsic colour redder than reality and too low 
a value of extinction. 

Field 6134 ((/,/?) = (81.46,-0.72)) is outside Cyg 0B2, but 
due to its relative proximity should exhibit similar properties to 
fields 5985 and 6010 at distances less than that to Cyg 0B2. Fig[23l 
demonstrates that the local (< 1 kpc) behaviour of extinction in 
fields 5985, 6010 and 6134 is much the same. This is as it should 
be, since at these near distances the field centres are less than 50 pc 
apart and likely to be sampling much the same stellar environment. 
Further examination of the results for field 5985 reveals a 30% in- 
crease in the proportion of objects in luminosity classes IV and III, 
relative to field 6134, suggesting that MEAD is indeed capable of 
dealing with the brighter absolute magnitudes of PMS objects. 



6 DISCUSSION 

We have presented an algorithm, MEAD, and shown how it can be 
applied to the r\ i' and narrow-band Ha photometry forming the 
IPHAS database to map out extinction in the Galactic disc at high 
angular and depth resolution. The algorithm is fully empirical in 
that it depends only on absolute magnitude calibrations (especially 
the Houk et al 1997 HIPPARCOS main sequence) and not on any 
Galactic model. Because of this it is better able to deal with per- 
turbations away from a smooth structure. This is crucial as it is 
these perturbations which are of real interest, in view of their pos- 
sible origin due to either accreted satellite galaxies or to the Milky 
Way's still-debated spiral structure. 
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Figure 25. The retrieved extinction-distance relationship for a 10' x 10' box 
in the direction {l,b) = (180,0) using the final version of MEAD (black 
dots), near MS early-A stars only (green open squares) and by assuming 
that all objects are dwarfs (red crosses). 



already presented a preliminary catalogue of emission line stars ex- 
tracted from above the main stellar locus. 

In the future, MEAD could also be applied to the results of the 
upcoming VST/OMEGACAM Photometric Hydrogen Alpha Sur- 
vey of the Southern Galactic Plane (VPHAS-i-). Combining the re- 
sults from VPHAS-i- and IPHAS would lead to an extinction map 
of the entire Galactic plane in unprecedented detail. VPHAS-i- and 
UVEX extend the wavelength coverage of IPHAS to the u' and g' 
bands over the full plane. These additi onal bands open up the inves- 
tigation of variations in metallicity fcf. llvezic et aLboOSh and in the 
extinction law. Here, in limiting ourselves to red wavelengths, we 
have been able to work with a standard R = 3.\ law as reasonably 
representative of most sightlines. Combination with near-infrared 
wavelengths is also becoming an option as the UKIDSS GPS pro- 
gresses and VISTA (Visible and Infrared Survey Telescope for As- 
tronomy) begins data collection. 

As a next step, the impending completion of IPHAS observing 
creates the opportunity to apply MEAD to a full 3D mapping of the 
northern Galactic Plane (Sale et al. in preparation). 



ljuric et aP (l2008h are able t o successfully c orrect for extinc- 
tion using the 2D maps due to 'Schlegel et al.! (Il998 ). since the 
majority of their observations sample high Galactic latitudes. For 
studies of the Galactic disc at low latitudes, a more sophisticated 
treatment of extinction is required. The algorith m MEAD pr e sented 
here is a means to achieving this. While the lJuric et iD (l2008h 
SDSS sample is effectively volume limited for most spectral types, 
leading to main sequence stars dominating the sample, the IPHAS 
database is more clearly magnitude limited, requiring the disentan- 
gling of luminosity class. This capability is a critical new feature of 
MEAD. 

The new power of MEAD can be illustrated by comparing the 
di stance-extinction relationships obtained (i) using A stars only (as 
in iDrew et aLll2008h . (ii) analysing intrinsic colour, whi lst classify- 
ing all stars as dwarfs, (iii) using MEAD in full (Fig. [25] ). lDrew et aP 
(l2008h demonstrated that extinctions and distances of near Main 
Sequence early-A stars were particularly straightforward to derive 
from their IPHAS colours and apparent magnitude. This approach, 
on its own, could in principle be used to map extinction, but near 
Main Sequence early A-stars represent only a small proportion of 
objects in the IPHAS catalogue. By accessing a broader range of in- 
trinsic stellar colours MEAD is able to map extinctions on finer spa- 
tial scales, with greater accuracy and over a wider distance range. 
This figure also demonstrates why it is so important to analyse lu- 
minosity class as well as intrinsic colour. 

As well as estimating the distance-extinction relation, MEAD 
also returns distance, intrinsic colour and extinction estimates for a 
large proportion of objects in the IPHAS database. Interpreting the 
resultant catalogue, will be non-trivial, but will enable the structure 
of the Galactic plane to be studied in a manner similar to the anal- 
ysis of the Milky Way at higher latitudes by Juric et al. ( 2008). For 
example, scale heights and lengths of both the dust and stars within 
the disc will be derivable. 

A side benefit of MEAD is that it generates lists of objects 
that do not conform to the most common spectral classifications. 
These 'rejected' objects may still be of interest for focused studies 
of rare phases of stellar evolution. Whilst, most of the objects in this 
category will be late K or M -type stars, or noise dominated, others 
would be extreme red giants I Wright et al.l2008l). sub-lumino us ob- 
jects fSection [33] ) or emission line stars. Witham et al.! (l2008l) have 
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